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Abstract Cosmologies including continuous matter creation are able to reproduce the 
main properties of the standard /1CDM model, in particular in cases where the particle 
and entropy production rates are equal. These specific models, characterized by a mass 
density equal to the critical value, behave like the standard /1CDM model at early times 
whereas their late evolution is similar to the steady-state cosmology. The maximum 
amplitude of density fluctuations in these models depends on the adopted creation 
rate, related here to the parameter fi v and this limitation could be a difficulty for 
the formation of galaxies and large-scale structure in this class of universe. Additional 
problems are related with predictions either of the random peculiar velocities of galaxies 
^-( or the present density of massive clusters of galaxies, both being largely overestimated 

I ■ with respect to observational data. 
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1 Introduction 

vo' 

\l ■ The standard model in cosmology, the so-called /1CDM model, assumes the presence 

,—jl. , of a constant cosmological term in Einstein's equations and that the universe is con- 

stituted, besides baryons, photons and neutrinos, by a dominant weakly interacting 
component of unknown nature dubbed dark matter. This model gives the best descrip- 
tion of the present data as, for instance, those resulting from the analysis of the seven 
years operation of the WMAP satellite [T] and supernova distances [2J. 

The inclusion of the cosmological term in Einstein's equations represents the most 
simple and natural explanation for the observed acceleration of the expansion of the 
universe detected through supernova data. However, arguments against the inclusion 
of such a cosmological term are often found in the literature. The first one is the 
so-called "coincidence" problem, which can be stated in the following way: why do 
we observe presently an almost equal contribution of matter and the cosmological 
term to the total energy budget of the universe, considering that these components 
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evolve differently with time? Such a "coincidence" would suggest that we are living 
in a particular moment of the history of the universe, contrary to a "cosmological 
principle" stating that we are not in a special place in the universe either in space or in 
time. Notice that if the accelerated phase started around z w 0.75, such a phase covers 
about 48% of the existence of our universe! The second argument is related to the 
interpretation of the cosmological constant as the vacuum energy density. As it is well 
known, formal quantization of classical fields leads to a divergent energy density for the 
vacuum state, which can be avoided by imposing a physical cut-off [3]. If we adopt as 
a cut-off the Planck scale, then the theory predicts an effective cosmological constant 
corresponding to an energy density that is about 120 orders of magnitude higher than 
the observed value. Convincing arguments against these difficulties were recently given 
by Bianchi and Rovelli [3] , who concluded that the "coincidence problem" is ill defined 
and that the identification of the cosmological constant with the vacuum energy density 
is probably a mistake. 

Face to this debate on the existence (or not) of the cosmological term, alterna- 
tive models have been discussed in the literature [5], IS], [7], El, [2]. In some of these 
models, the acceleration of the expansion is driven by a negative pressure term asso- 
ciated to particle creation at the expense of the gravitational field, an original idea 
proposed by Zeldovich [10] almost 40 years ago. Particular cases in which the particle 
production rate is proportional to the Hubble parameter give results similar to the 
canonical ylCDM model [TT|, [T2]. From a thermodynamic point of view, particle pro- 
duction from the gravitational field can be considered in the context of open systems, 
where the "heat" received by the system is due entirely to the change in the number 
of particles [13]. In the standard model, matter is produced suddendly at the initial 
singularity and during the re-heating at the end of the inflationary phase. After that 
the universe evolves adiabatically, i.e., both the entropy per unit of comoving volume 
and the entropy per particle remain constant during the expansion. 

In the present paper we review the properties of specific models including particle 
creation, emphasizing some particular aspects not previously noticed. We investigate 
also the linear growth of density perturbations in this class of cosmological models and 
we will show that, depending on the creation rate, the amplitude of the density contrast 
is unable to grow beyond a certain limit. The critical creation rate below which density 
fluctuations are able to grow continuously up to present time is marginally compatible 
with the parameters of the ACDM moldel. This could represent a difficulty for the 
formation of galaxies and the large-scale structure, condemning this particular class 
of cosmology. This paper is organized as follows: in Section 2 the principal physical 
aspects of the model are revisited; in Section 3 the growth of density perturbations is 
examined and finally, in Section 4 the main conclusions are given. 



2 The cosmological model 

In the present investigation, for simplicity, we will neglect the evolution of the baryonic 
component and we will focus our attention on the dominant (dark matter) component 
only. The number of baryons per unit of comoving volume remains conserved during the 
expansion, i.e., n^a = constant, where nj, is the comoving baryon number density and 
a — a(t) is the scale factor. On the other hand, the bulk of the dominant component is 
produced as in the standard scenario but now the possibility of a continuous creation 
of dark matter particles is included. This implies that we renounce to the common 



idea that the expansion of the universe could be described by "closed" system and we 
assume, following Prigogine et al. [13] that, in fact, the expansion is described by an 
"open" thermodynamic system. Under these conditions, the first law can be written as 

TdS = dE + PdV - (jdN (1) 

where the chemical potential p is here associated to the non conservation of the parti- 
cle number. Introducing respectively the entropy density s, the energy density p, the 
particle number density n, the enthalpy density h = (P + p), replacing into eq. [TJ 
developing and using the Euler's relation 

np = (P + p) - Ts = h - Ts (2) 

one obtains after some algebra 

fdp hdn\ = / digs dlgn \ .. 

V dt n dt ) S \ dt dt ) ( ' 

which is essentially the Gibbs relation for the system (see [TJ] for a covariant derivation 

of this relation). 

Solutions of eq. [3] permit the study of different particle production scenarios but 

relations either between the state variables s and n or for their production rates are 

required, since they define the model. The simplest possibility, considered either by |11| 

or []2] , corresponds to the case where the entropy per particle remains constant during 

the expansion process, i.e., s/n = constant. This hypothesis implies that the energy 

density a function of the particle number density only and not of the temperature 

also. Moreover, it is worth mentioning that such a condition is not equivalent to an 

adiabatic expansion as it occurs in the standard model since, as we shall see later, there 

is entropy production and the aforementioned condition expresses only the fact that 

the relative rates of particle and entropy production are equal. In this case, the right 

side of eq. [3] is zero and the particle production rate is related to the energy density 

variation rate by 

dp h dn 

IE ~ nUt ^ ' 

Let the stress-energy tensor of the dark matter fluid be 

Tab = (p + n)u a u b - n 9ab (5) 

where 77 = P(p) + P c is the effective pressure acting on the fluid, with the first term on 
the right side representing the pressure due to kinetic motions and interactions between 
particles and the second, P c , being a new term associated to the particle production. 
Considering a Friedman-Robertson- Walker (FRW) metric for a spatially flat spacetime, 
i.e., 

ds 2 = -dt 2 + a 2 {t) {dr 2 + r 2 dn 2 ) (6) 

from the equation T* . fe = 0, one obtains for the variation rate of the energy density 

g + 3f(* + ft) = (7) 

Then, from eqs.[4]and[7] one obtains immediately an expression for the pressure asso- 
ciated to particle production 



where the Hubble parameter H was defined as usually, namely, H — a/a. It should be 
emphasized again that the above relation is valid only if the ratio s/n remains constant 
in the process. As expected, the relation above indicates that if the particle produc- 
tion rate is positive in a expanding universe (H > 0), the pressure P c is negative, 
contributing to accelerate the expansion. The simplest world model without a cosmo- 
logical constant, but including a negative pressure term due to particle production is 
obtained from the ansatz P c = —A, where A is a positive constant having the dimen- 
sion of an energy density. This ansatz is equivalent, as we shall see below, to the usual 
assumption that the particle production rate is proportional to the Hubble parameter 
if the created particles are "cold" and interact "weakly". Consequently, the enthalpy 
density is simply h = p and the relation between the particle and the energy densities 
is given by p = nrnc , where m is the rest mass of the created particles (supposed to 
be uncharged). Under these conditions, eq. [8] can be recast as 

dn „ „ 3A „ ,„. 

— +3Hn= rrH (9) 

at mc z 

Notice that now the particle conservation equation has a (positive) source term pro- 
portional to the Hubble parameter as mentioned above. Using the scale parameter a 
as the independent variable instead of time, the equation above can be rewritten as 

dn „ 3A ,_„. 

+ 3n = j (10) 



dig a mc 

and similarly for the energy density 

dp 



+ 3p = 3A (11) 

dig a 

Integration of eq. [TT]is trivial and is given by 

p = A+ (^l (12) 

where p* is the energy density when a = 1 (the present time). The equation above can 
be recast as 

p = p* (ftv + Q m a~ ) (13) 

which is formally identical to the expression of the energy density in the case of the 
standard /1CDM model, if one identifies O v = A/p* and Q m = 1 — Q v . However, in the 
present case fl v is not the density parameter associated to the cosmological term but 
to the creation rate. Similarly, O m in the present case is not equivalent to the matter 
parameter density as in the standard model, since here the present matter density 
is equal to the critical density. Integration of eq. \TU\ is also trivial and gives for the 
evolution of the particle number density 

n = n* [Q v + n m a~ 3 ) (14) 

where n* is the present particle number density. For high redshifts or a << 1, the 
particle number density satisfies the condition na = n*Q m — constant, similar to the 
standard model. However, in the future, when a » 1, the standard model predicts that 
either for baryons or for dark matter, the particle number density goes to zero due to 
the expansion of the universe while the present model predicts a constant density equal 
to n = n*O v for dark matter particles. In other words, the present model predicts a 



future state of the universe constituted by dark matter only since the density ratio n^/n 
goes to zero. This future behaviour of dark matter is identical to that of the steady- 
state cosmology, since when a » 1 the creation rate will be equal to the expansion 
rate. It worth mentioning that the present matter density in this model is higher than 
the standard ylCDM by a factor of l/fi m , a fact with observable consequences, as we 
shall see later. In reality, these are not the only differences with the standard model. 
Since the ratio s/n remains constant during the expansion, the evolution of the entropy 
density is given by 

s — s* (j?!, + Q m a~ ) (15) 

where s* is the present entropy density of dark matter. The equation above says that the 
evolution of the entropy density also differs from the standard model, with a temporal 
behaviour similar to that of the particle density or, in other words, in the future the 
entropy of the universe will essentially dominated by that of dark matter. Notice also 
that in the standard model, the entropy in a proper volume corresponding to a unit 
comoving volume (sa ) is constant but in the present model there is entropy production 
at a rate 

^l=3 St ,n v a 3 H (16) 

Once the past behaviour of the present model is similar to that of the standard model, 
both models satisfy the different tests (supernova distances, baryon acoustic peak, 
"shift" parameter associated to the CMB, age of the universe) if the parameters of 
both models are conveniently interpreted (see, for instance, reference 12). 

The particular cosmology here considered has a deceleration parameter q formally 
similar to the standard model, i.e., 

_ aa _ 1 (f2 m a~ 3 — 2Q V ) 
q ~ a 2 ~ 2 (j?„ + n rn a- 3 ) [ ' 

Thus, the present value of this parameter is simply 

g* = ~(l-3i?„) (18) 

Since observations of supernova distances indicate O v = 0.7 ± 0.1 ([IS]), from the 
equation above it results that g* = —0.55 ± 0.15. 



3 The evolution of density perturbations 

In the present study, in order to obtain the linearized equations, the Newtonian ap- 
proximation will be used. This approximation is justified in the weak field limit, i.e., 
when the velocity of peculiar motions satisfies the condition V p « c and the scale of 
the perturbations is much less than the Hubble radius (see ref. [H]). Notice that in the 
case of the "standard" model, a fully relativistic treatment leads essentially to same 
results obtained through the weak field approximation. 

The first equation is the particle number conservation (eq.[9]), which can be recast 
as 

^ + V • (nU) = aH (19) 

where a = 3A/mc = ?>Q v n*. The velocity field U includes the Hubble flow and pe- 
culiar velocities resulting from gravitational interactions between density peaks. The 



second equation expresses the gravitational aceleration in terms of the mass distribu- 
tion, i.e., 

V • g = 4nGmn (20) 

Finally, the equation of motion, which in absence of pressure gradient terms can be 
simply written as 

f + u.f +g = » 

These equations are similar to those in the usual Newtonian approximation except by 
the fact that now a source term is present in the continuity equation. Defining the 
perturbed quantities n\, gi and V p by the relations, 

n(r,i) = n (t) + ni(r,t) (22) 

g(r,t) = go(r,t)+ gl (r,t) (23) 

U(r,t) = #r + V p (r,t) (24) 

introducing the density contrast S — m/no, replacing the perturbed quantities into 
eqs. 1191 l20l 1211 expanding and retaining only first order terms, one obtains respectively 
for the linearized continuity, force and motion equations 

f «-~ -.t 4nGm . .„_. 

8 + V-V p = — — a5 (25) 

3 H 

V ■ gi = 4itGmn 8 (26) 

and 

-gf- + H V P + gl = (27) 

It should be emphasized that here the subscript "0" denotes unperturbed values while 

quantities taken at present time (z — 0) are denoted by the subscript "* " . Writing the 

perturbed quantities in terms of the Fourier components, the system of linear equations 

above becomes 

x ik ■ V fc AirGma 

Ok H = TFT — °k (28) 

a oriQ 

tkgfc = 4nGmn S k (29) 

-gf+H o V k + Sk =0 (30) 

After some lengthy but trivial algebra with these equations, one obtains finally for the 
evolution of the density contrast in the linear approximation 

4 + (2^0 + ^|p) h + 47rGm [a (l - |±) - n ] 5 k = (31) 

Notice if there is no particle creation (a = 0) the usual equation describing the evolution 
of the density contrast in absence of pressure gradients is recovered (see, for instance 

El). 

In order to solve numerically eq. 1311 it is convenient to introduce a dimensionless 
time variable r = t/t*, where t* is the present age of the universe. Using now the 
notation d5 k /dr = 8 k and d S k /dr = S k , eq. [31] can be written in a general form 

h(j) + 9l(r)S k (T) + g 2 (r)8 k (r) = (32) 



where explicit expressions for the functions g\ (r) and 32 (t) depend on the adopted cos- 
mological model. For the standard yICDM model, these functions are given respectively 
by 

Si (r) = 2/3* f(a(r)) = 2/3* ^JT+^FHj) (33) 



and 



S2(r) = -|/3» 2 X a- 3 (r) 



(34) 



In these equations, /3* = H*t*y/f2 v , \ = Q m jfi v and the scale factor as a function of 
the dimensionless time is 



o(t) = x 1/3 5m/i 2/3 (3/3»r/2) 



(35) 



When particle creation is included (a / 0), the functions <?i(r) and 52 (t) become 



ff i(r) = 2/3*/(a(r)) 



1 + 



4/2(o(t)) 



and 



52 (t) 



I* 



(4 + 7 X a' 3 (r)) 2 



2(l+ X a- 3 (r)) 



/Vm) 



(36) 



(37) 



In both cases the parameter /3* can be estimated from eq. 1351 using the condition 
o(l) = 1. Thus, /3* = 2ArcSinh(x~ )/3, implying that the only free parameter in 
these equations is Q v - 

Numerical solutions of these equations were obtained by adopting initial conditions 
appropriate to the decoupling between matter and radiation at z ~ 1100 and using 
the fact that at these high redshifts both models are similar to the Einstein-de Sitter 
cosmology. Moreover, since there is no continuous production of baryons, the transition 
from an "opaque" to a "transparent" universe should occur at the same redshift as in 
the canonical model. Thus, the evolution of the density contrast when a « 1 is 
8(t) oc a(t). 

Our results indicate that for models with a/0, i.e., including particle creation, 
the density contrast increases, reaches a maximum before the present time and then 
decreases if the parameter il v is higher than a critical value equal to f2 v . CT it = 0.666. 
For this particular value, the maximum occurs at the present time and for lower values, 
the maximum occurs in the future, i.e., when a > 1. This behaviour is illustrated in fig. 
1 where the evolution of the density contrast as a function of the scale factor and for 
different values of the parameter Q v is shown. For a < 0.03 orz> 35, as expected, all 
models coincide since they follow the same behaviour of the Einstein-de Sitter solution. 

The existence of a maximum in the evolution of the density contrast is a conse- 
quence of a fundamental difference between models including particle creation and the 
standard ylCDM cosmology. These differences are clearly seen in a close inspection 
of the coefficients appearing in eq. [31] either when a = (vlCDM model) or when 
a/0, particle creation model. The first coefficient, corresponding to the second term 
on the left (or the coefficient of the first derivative of the density contrast) is simply 
equal to 2H — 2a/ a in the case of the standard model. This is a damping mecha- 
nism related to the expansion of the universe. When particle creation is included, this 
damping coefficient increases from 2H up to 2.36H for a = 1 and Q v = 0.73. Such 
a variation is a small effect not responsible for the existence of a maximum, although 



Models with matter creation 




0.2 0.3 0.4 0.5 0.6 

scale factor a(t) 



Fig. 1 Density contrast evolution as a function of the scale parameter. Labels indicate different 
values of the parameter fi v 



it contributes to decrease slightly the amplitude of the density contrast. The second 
coefficient, corresponding to the third term on the left of eq. 1311 can be recast as 



47t Gm 



'i ( 1 + - ) - ih) 



(38) 



where q is the deceleration parameter. In the standard model (a — 0) the coefficient 
above is negative, representing the well known gravitational instability. When parti- 
cle production is included, an extra term counterbalancing gravity appears and the 
amplitude of the density contrast increases as far as the condition 

2 < 1 + * a ~ 3 > a > 1 (39) 

(4 + 7 X a- 3 ) 

is satisfied which, for adequate values of x — dm/^v, is equivalent to the condition 
a < 0.5. Consequently, depending on the value of H v , a maximum amplitude is attained 
in the late evolution of the density contrast, whose value depends on the particle 
creation rate. 

In fig. 2, the density contrast evolution for Q v — 0.73 is shown for both the standard 
and the particle creation model. The density contrast of the later at the present time 
is about 30% lower than that derived for the standard model, with the maximum 
occurring at z ~ 0.2. 

The effective damping of the density contrast induced by particle production poses 
some difficulties to the process of galaxy and large-scale structure formation by grav- 
itational instability. However, other problems related with the growth of the density 
contrast exist in this model. For instance, the present root mean square value of the 
peculiar velocity of galaxies in the linear theory when particle creation is absent is 
given by [IT] 

<Vp 2 > 1 ^Ha(^-)^T 2 (40) 

V dl 9 a J 




Fig. 2 Comparison between the density contrast evolution resulting from the standard TlCDM 
model and the present cosmology including matter creation. 



where J2 = 150h~ Mpc is the second moment of the two-point correlation function of 
galaxies [18] . The introduction of the particle creation process modifies the continuity 
equation, as we have seen previously, since a source term is now present. Consequently, 
the resulting equation for the root mean square velocity of galaxies in the linear theory 
is now given by 



< V p 2 > 1 ' 2 -- 



Hf(a)a 



1+^ 



ill) 



(41) 



2 n f(a) 

where we have defined as in reference [T7] f(a) — dig 5 /dig a. Numerical solutions 
of the equation above indicate that the present root mean square velocity practically 
independs on the parameter Q v , being equal to ~ 1200km/ s. This is considerably 
higher than the value derived from observational data, i.e., 325fcm/s |19| . Moreover, 
as mentioned before, the present matter density is greater than that of the standard 
ylCDM model by a factor of Qm . As a consequence, the predicted density of clusters of 
galaxies above a given mass is higher either than predictions of the standard model or 
observational data. For clusters with masses higher than 10 Mq, data by Bahcall and 
Cen 20 indicate that the density is about 7.2 x 10 Mpc while from the creation 
model a density of about 3.8 x 10 Mpc is predicted. 



4 Conclusions 



Cosmological models described by "open" thermodynamic systems, i.e., including par- 
ticle creation at the expense the gravitation field can reproduce formally the past 
evolution of matter, mimicking the presence of a cosmological constant in Friedman 
equations. This class of models results from the assumption that entropy and particle 
production rates are equal and that the negative pressure resulting from the particle 
creation process is a constant. This last hypothesis is equivalent to the assumption of 
a production rate proportional to the Hubble constant. It is worth mentioning that 
in the steady-state cosmology the creation rate required to compensate the expansion 
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and to maintain a constant particle density is r ss = 3H*n*, which should be compared 
with the rate of the present model, given by r cre = ?>Q v H*n*. Thus, the present model 
has two asymptotic behaviours: in earlier times behaves as the standard model or the 
Einstein-de Sitter solution and in late times behaves as the steady-state cosmology, 
except by the fact that the baryon-to-dark matter particle density ratio goes to zero. 
Difficulties appear for the creation model when the evolution of density fluctuations 
is regarded closely. In this model, the particle production mechanism mimics also a 
damping mechanism that limits the growth of the density contrast, constraining the 
process of galaxy and large-scale structure formation. Moreover, the predicted root 
mean square peculiar velocity of galaxies in the linear theory is about four times higher 
than observations. The creation model has also difficulties to explain the observed 
density of massive clusters of galaxies, predicting presently (z — 0) an excess of objects 
by a factor of fifty with respect to observations, consequence of the fact that the 
present matter density in this cosmology is equal to the critical value. Nevertheless 
it is important to mention that the present analysis is based on a linear approach 
and that the inclusion of non-linear terms, affecting the late evolution of the density 
contrast may eventually change our conclusions. Cosmological simulations including 
particle creation will be the next step to investigate these non-linear effects, that will 
be able to corroborate or not the present study. 
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